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ABSTRACT 


This  paper  considers  the  lot-sizing  problem  in  a multistage 
inventory  system.  External  demand  may  occur  at  any  stage,  and  is 
assumed  to  be  known  over  a finite  horizon.  A heuristic  iterative 
procedure  is  proposed  and  tested  for  finding  a periodic  review 
schedule  to  minimize  inventory  and  setup  costs. 
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1.  Introduction 

In  this  paper  we  consider  a dynamic  lot-sizing  problem  for  a multi- 
stage inventory  system.  The  lot-sizing  problem  is  to  determine  replenish- 
ment quantities  for  an  inventory  system  so  as  to  satisfy  all  demand  require- 
ments at  minimum  system  cost.  A multistage  Inventory  system  is  a connected 
set  of  stages  representing  the  steps  for  assembly  and/or  distribution  for 
a series  of  products.  Typically  we  may  characterize  a multistage  inventory 
system  as  an  acyclic  network  modified  to  include  links  between  the  stages 
and  the  outside  customer  demand;  examples  of  such  systems  are  given  in 
Figure  1.  A standard  classification  is  to  distinguish  between  assembly 
and  arborescent  inventory  systems;  for  an  assembly  system  (Figure  lb),  each 
stage  either  has  a unique  immediate  successor  or  satisfies  directly  outside 
customer  demand,  while  for  an  arborescent  system  (Figure  lc) , each  stage 
has  either  a unique  immediate  predecessor  or  no  predecessors  at  all.  A 
serial  system  (Figure  la)  is  both  an  arborescent  and  an  assembly  system, 
while  a general  multistage  system  need  not  be  either  an  arborescent  or  an 
assembly  system. 

The  objective  of  this  paper  is  to  propose  and  test  a heuristic  lot- 
sizing procedure  for  a general  multistage  discrete-time  inventory  system. 

We  assume  that  outside  customer  demand  is  known  by  period  over  a finite 
horizon,  and  must  be  satisfied  from  on-hand  inventory  with  backordering  not 
allowed.  All  lead  times  for  supplying  a stage  from  an  immediate  predecessor 
are  given,  and  without  loss  of  generality  are  assumed  to  be  zero.  A lot- 
sizing  procedure  specifies  a periodic  schedule  indicating  the  size  and 
timing  of  inventory  replenishments  so  that  each  stage  satisfies  the  demand 
placed  upon  it  by  succeeding  stages.  The  optimality  criterion  is  to  mini- 
mize total  cost  over  the  finite  horizon  where  there  are  two  types  of  costs 
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1:  Examples  of  Multistage  Inventory  Systems 
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at  each  stage:  a fixed  cost  for  placing  an  order,  and  an  inventory  holding 
cost  which  we  assume  to  be  proportional  to  the  end-of-period  inventory  at 
each  stage. 

There  has  been  a great  amount  of  work  done  on  dynamic  lot-sizing 
problems,  beginning  with  the  classical  work  of  Wagner  and  Whltln  [15]  for 
the  single-stage  problem.  Significant  extensions  have  been  made 
by  Wagner  [14]  to  Include  dynamic  cost  functions,  by  Zabel  [16],  Eppen , 
et.  al.  [5],  Kunreuther  and  Morton  [7], [8],  Blackburn  and  Kunreuther  [2],  and 
Lundln  and  Morton  [10]  to  extend  and  generalize  the  planning  horizon  theorems, 
and  by  Zangwlll  [17]  to  Include  backlogging  of  demand. 

For  the  dynamic  multistage  lot-sizing  problem,  Zangwill  [17]  has  shown 
that  for  a series  system  the  optimal  solution  is  contained  in  the  set  of 
extreme  flows  of  a single-source  network,  and  has  given  a dynamic  programming 
recursion  for  calculating  the  optimal  policies.  Love  [9]  also  considers  a 
series  Inventory  problem,  and  gives  an  alternate  algorithm,  which  exploits 
the  nested  property  of  the  solution. 

Veinott  [13],  Crowston  and  Wagner  [4],  and  Kalymon  [6]  have  considered 
more  general  multistage  problems,  and  have  given  algorithms  for  solving 
specific  problems;  however  in  all  three  papers,  the  algorithms  are  quite 
complex,  with  the  amount  of  work  required  increasing  exponentially  with  the 
number  of  stages,  or  with  the  number  of  time  periods,  or  with  both. 

Due  to  the  complexity  of  the  general  multistage  problem,  many  heuristic 
procedures  have  been  proposed.  The  most  common  form  of  heuristic  is  to 
consider  the  stages  sequentially,  starting  with  the  lowest  echelon  stages1, 
and  scheduling  each  stage  with  a single-stage  procedure  which  may  by  Itself 
be  a heuristic.  Examples  of  such  heuristics  are  given  by  McLaren  [11], 


1 


The  lowest  echelon  stages  are  those  nearest  to  the  customer. 


McLaren  and  Whybark  [12],  Biggs,  et.  al.  [1],  and  Blackburn  and  Millen  [3]. 

For  an  N-stage  system,  the  amount  of  work  necessary  for  these  heuristics, 
which  we  term  as  single-pass  heuristics,  is  comparable  to  that  needed  for 
solving  N single-stage  problems.  All  of  the  reported  work  has  been  restricted 
to  assembly  systems.  From  this  work,  a heuristic  given  in  [11],  [12]  seems 
to  give  the  best  performance,  although  any  conclusions  are  limited  by  the 
scope  of  the  computational  studies;  this  heuristic  is  discussed  in  greater 
detail  in  Section  3. 

The  intent  of  this  paper  is  to  present  and  test  a new  type  of  heuristic, 
a multipass  heuristic  for  the  multistage  lot-sizing  problem.  Whereas  a 
single-pass  heuristic  sequentially  schedules  each  stage  and  then  stops,  a 
multipass  heuristic  does  not  stop  after  the  "single-pass",  but  continues  to 
revise  the  current  schedule  in  an  iterative  fashion  until  no  further  improve- 
ments in  the  schedule' are  possible.  Admittedly,  the  multipass  heuristics 
require  more  computational  effort  than  single-pass  heuristics,  but  hopefully 
much  less  than  an  optimal  algorithm.  The  schedule  performance  of  a multi- 
pass heuristic  should  likewise  be  bounded  by  that  of  a single-pass  heuristic 
and  an  optimal  algorithm. 

The  remainder  of  the  paper  is  organized  as  follows:  In  the  next  section, 
a multipass  heuristic  is  developed  and  is  shown  to  be  monotonic  and  convergent 
In  Section  3,  we  report  and  discuss  our  computational  experiments  comparing 
the  multipass  heuristic  with  two  single-pass  heuristics,  and  with  an  optimal 
algorithm. 
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2.  The  Multipass  Heuristic 

The  multipass  heuristic  consists  of  two  phases,  one  in  which  the  current 
schedule  is  revised  and  the  other  in  which  various  stages  are  "collapsed" 
into  other  stages  based  on  the  current  schedule.  The  presentation  in  this 
section  will  first  develop  the  logic  for  revision  and  then  will  consider  the 
collapsing  procedure;  this  development  is  motivated  first  by  considering 
simple  systems  and  then  is  generalized  for  more  complex  systems. 


2.1  A Revision  Procedure  for  a Two-Stage  System 

Consider  a two-stage  system  as  depicted  in  Figure  2.  Item  2 is  the  single 
component  for  item  1,  the  final  product.  For  each  unit  of  item  1 assembled, 

6 units  of  item  2 are  required.  This  two-stage  structure  is  the  simplest 
multistage  system,  and  is  used  here  in  order  to  simplify  the  presentation 
of  the  proposed  heuristic.  All  results  for  this  system  can  be  generalized 
to  more  complex  multistage  structures. 


©-S-© 


customer  demand 


Figure  2:  A Two-Stage  System 


Define  the  following: 

d^t  « external  demand  for  item  i in  time  period  t, 

S^  - setup  cost  for  item  i, 
h^  » holding  cost  for  item  i, 

Vjt  “ variable  production  cost  for  item  i in  time  period  t, 

T ■ length  of  scheduling  horizon  in  periods. 

We  assume  that  the  initial  inventory  for  both  items  is  zero,  the  production 
lead-time  is  zero,  and  d^  > 0.  Also,  note  that  item  2 may  have  external 


demand,  independent  of  the  demand  imposed  on  it  by  item  1.  For  notational 
convenience,  we  will  identify  the  single-stage  problem  with  demand  require- 
ments dt  for  t-l,...,T,  setup  cost  S,  holding  cost  h,  and  variable  produc- 
tion cost  vt  for  t“l,...,T,  as  VJWld^.S.h.v^ , i.e.  the  Wagner-Whitln  problem. 
Now  the  initial  heuristic  (HI)  is  as  follows: 

(HI): 

1.  Solve  WW[d^t ,S^,h^,v^t ] . If  x^t  is  the  production  quantity  of 
item  1 in  time  period  t,  then  B is  the  demand  imposed  on  item  2 
by  the  production  schedule  of  item  1. 

2.  Set  * d^j.  + “ total  demand  for  item  2 in  period  t.  Solve 

WW[d2t,S2,h2,V2t].  From  this  solution,  determine  m marginal 
cost  of  increasing  demand  for  item  2 in  period  t by  one  unit. 

3.  Set  Vjt  “ v^t  + BY2t  = restated  variable  production  cost  in  period  t. 
Solve  WW[d^t ,S^,h^,v^t] . If  there  is  no  change  in  the  schedule  for 
item  1,  then  stop;  if  there  is  a change,  then  return  to  step  2. 

This  procedure  is  an  Iterative  procedure  in  which  two  single-stage 
problems  are  solved  at  each  iteration.  We  will  show  that  the  procedure  is 
nonitv-reasing  in  cost,  and  hence,  convergent.  First,  though,  we  indicate 
how  the  marginal  costs  Yj,-  are  determined. 

Define  Y2t  to  marginal  cost  of  increasing  by  one  unit  the 

demand  requirements  of  item  2 in  time  period  t given  a production  schedule 
for  item  2;  Y2t  can  he  interpreted  as  a shadow  price  on  demand  for  a fixed 
schedule  for  item  2.  Let  T be  the  last  time  period  prior  to  t in  which 
there  is  production  of  item  2 (i.e.  the  setup  cost  is  incurred).  That  is, 
if  X2t  is  production  in  time  period  t,  then 
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x2t  > 0 

X2J  ' 0 


for  j-T+1, . . . ,t , 


or  x^t  > 0 and  t “ T. 


Now  we  specify  Y-  as 

i.  C 


v2t  + (t-x)"h2 


That  is,  if  demand  in  period  t is  increased  by  one  unit,  production  in 
period  T (x2t)  increases  by  one  unit  with  an  incremental  cost  consisting  of 
the  variable  production  cost  in  period  t,  plus  the  holding  cost  from 

period  T to  t which  is  (t-T)*h.  Note  that  we  assume  x always  exists;  this 
is  true  provided  that  there  is  no  initial  inventory  and  demand  in  the  first 
period  (external  or  induced)  is  nonzero. 

We  now  show  that  the  procedure  coverges.  Define  the  following: 
cln  ■ cost  (using  original  costs  S^.h^.v^)  for  schedule 
of  item  i after  the  n1"*1  iteration, 

{xj  } ■ schedule  for  item  i after  nth  iteration, 

{y2t}  ■ marginal  costs  generated  after  n1*1  iteration  from  {*2tK 
Axt  - x]*^  ~ “ change  in  schedule  for  item  1 in  period  t. 

At  the  n+1  iteration,  for  item  1 we  solve  the  "revised"  problem 

,-.r  . „ . /v  , , * . n n . ^ ...  r n+1, 

WW t^it*si»^i«vit where  v^  ” v^t  + 8Y2t»  to  obtain  the  schedule  (x^t  i. 

If  the  cost  for  schedule  {x"**}  using  the  actual  production  costs  v^t  is 

c^  n+^,  then  the  "revised"  cost  using  v^t  is  c^  + Z y2t (fix"**) • 

Since  {x"^}  is  the  optimal  schedule  for  the  "revised"  problem 


WW[dlt,S1,h1,vlt]  , we  have 


'l.fl  * j,  ^<'1  i cl„  + j,  *?,«■$,> 


t 
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where  the  RHS  is  the  cost  for  the  "revised"  problem  using  the  schedule 
{*"}.  We  can  rewrite  (1)  as 

(2>  Cl,n+1  i Cl,n  - ^t(Btet> 

St 

For  item  2 at  the  n+1  iteration,  the  demand  is  given  by 
d"^  “ d2t  + 8x"^.  The  previous  schedule  {x2t}  ma^  '5e  m°dlfied  to  satisfy 
this  demand  {3”^}  by  just  adjusting  the  positive  order  quantities  (i.e., 
x2(.  >0);  if  {y  } is  the  modified  schedule,  we  have 


If  x2t  - 0, 


'■  ' X“t  + A( 


qn+l  ~r 
d2j  - d2 


x.  +8  Z Ax 


if  x‘‘t  > 0, 


where  T is  defined  such  that  x"j  = 0 for  j=t+l , . . . ,T-1 , and  x^T  > 0,  or 

T = T+l.  The  cost  for  item  2 for  the  schedule  {y  } can  be  shown  to  be  equal 

T C 

to  c,  +8  2 y”  (Ax  1 ; for  c„  being  the  cost  for  the  optimal  schedule, 

zn  ^t  t z,n+i 

we  have  by  definition 


c2,n+l  i c2n  + 6 * ^t(Axt>* 

t=l 


Adding  (2)  and  (3),  we  obtain 


(4)  c,  + c,  < c,  + c„ 

l,n+l  2, n+1  — l,n  2,n 

Thus  at  each  iteration,  total  cost  is  nonincreasing.  If  we  ignore  the  possi- 
bility of  cycling,  the  procedure  converges  since  the  single-stage  (WW)  problem 
considers  only  extreme  point  schedules,  of  which  there  are  a finite  number. 
Nevertheless,  it  may  be  possible  for  the  solution  to  cycle;  that  is,  the 
procedure  cycles  over  a finite  number  of  solutions  with  no  cost  improvement. 
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This  may  be  avoided  by  stopping  the  procedure  if  any  solution  is  repeated, 
or  by  perturbing  the  data  so  that  no  two  solutions  may  have  the  same  cost. 
One  possible  perturbation  is  to  redefine  the  setup  cost  for  stage  1 in  time 
period  t as  *=  + e1  where  c is  a small  positive  constant. 


2.2  Extension  to  More  Complex  Structures 

The  heuristic  (HI)  and  its  convergence  properties  are  extendable  to 
more  complex  structures.  For  instance  consider  the  following  two-echelon, 
6-stage  system  (Figure  3) : 


customer  demand 


customer  demand 


Figure  3:  Six-Stage,  Two-Echelon  System 


Define  0^  to  be  the  number  of  units  of  item  i needed  for  one  unit  of 


item  j.  Then  (HI)  is  modified  as  follows: 

1.  In  step  1,  the  independent  single-stage  (WW)  problems  for  the  first 
echelon  (items  1 and  2)  are  solved;  these  solutions  determine  the 
demand  for  echelon  II  items. 

2.  In  step  2,  the  independent  WW  problems  for  the  second  echelon 
(items  3,  4,  5,  and  6)  are  solved.  Note  that  for  item  5,  we  have 
d 


5t 


d^t  + B^x^t  + ®52x2t  w^ere  ^5ixit  *8  t*ie  *nduced  demand  from 


the  schedule  for  item  i,  i=l,2.  For  each  item  in  echelon  II,  the 


marginal  costs  y can  be  computed  as  before. 
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3.  Step  3 la  Identical  to  step  1 except  that  the  "revised"  WW  problems 
are  solved  for  Items  1 and  2.  Here  the  variable  production  costs 
are  adjusted  to  reflect  the  effect  of  the  higher  echelon.  For 
instance,  we  have  - vu  + + 6^  + 65^5,.. 

The  proof  of  the  convergence  of  this  procedure  is  identical  in  structure  to 
that  given  in  the  previous  section. 

The  extension  of  the  heuristic  to  more  than  two  echelons  is  slightly 
more  complex.  The  original  heuristic  is  inherently  geared  for  iterating 
between  two  echelons.  With  more  than  two  echelons,  we  must  decide  in  what 
order  the  procedure  cycles  over  the  echelons.  This  can  be  best  illustrated 
by  considering  the  following  three-echelon,  three-stage  system  (Figure  4): 


(D- 


Q- 


customer  demand 


Figure  4:  Three-Stage  Serial  System 


Here  each  stage  corresponds  to  an  echelon. 

Again  the  revisions  to  (HI)  are  straight-forward  for  the  structure. 
Demand  at  the  higher  echelons,  given  schedules  for  the  lower  echelons,  is 


2t 

d2t  + 

B21Xlt 

3t 

d3t  + 

632X2t 

The  revised  variable 


costs  are  computed  as 


where 

and 


It 


Vlt  + 


621Y2t 


Y2t  " °2t  + (t_T)h2 
C2t  “ V2T+e32Y3T 
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for  T being  the  last  production  period  of  item  2 prior  to  period  t.  Note 
that  l2t’  t*'e  mar8inal  cost  for  item  2,  reflects  not  only  the  incremental 
cost  at  stage  2,  but  also  the  Incremental  cost  at  stage  3. 

Now,  for  this  system  two  procedures  suggest  themselves  for  iterating 
over  the  echelons: 

a)  echelon  I -*•  II  III  -*■  I -*•  II  ... 

Here  a completely  cyclic  procedure  is  used;  having  scheduled  the 
echelons  from  I to  111,  the  procedure  returns  to  echelon  I and  tries  to 
improve  the  schedule. 

b)  echelon  I -►  II  •+  III  -*•  II  -*•  III  ...  (until  II  and  III  converge  ) -*•  I 
-*•  II  -*■  III  II  -►  III  ... 

Here,  given  a schedule  for  echelon  I,  the  two  echelon  problem  consist- 
ing of  II  and  III  is  completely  solved.  Given  the  solution  to  II  and 
III,  the  procedure  returns  to  I and  revises  its  schedule.  The 
procedure  then  repeats  with  the  new  schedule  for  I. 

It  is  not  clear  which  of  these  two  procedures  is  more  effective.  Preliminary 
computational  tests  on  small  examples  suggest  that  the  schedules  are  relative- 
ly insensitive  to  the  procedure,  but  that  procedure  b)  requires  more  work 
than  a) . 

Both  procedures  a)  and  b)  can  be  shown  to  converge.  The  convergence 
proof  for  procedure  a)  is  given  in  the  Appendix.  For  procedure  b)  the  proof 
follows  directly  that  given  for  the  two-stage  system;  to  see  this,  note  that 
this  procedure  may  be  viewed  as  the  solution  to  a nested  series  of  two-stage 
problems,  each  of  which  is  well-behaved. 

The  M-echelon,  N-stage  problem  (M  < N)  is  now  Just  a combination  of 
the  previous  two  extensions.  Again,  for  M > 2,  we  have  the  problem  of  deciding 


how  to  iterate  over  the  echelons. 
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2.3  Additional  Improvements  by  Stage  Collapsing 

The  heuristic  (HI)  Is  an  Improvement  heuristic.  Unfortunately  examples 
can  be  easily  constructed  to  show  that  It  Is  not  an  optimal  procedure. 

Indeed,  It  Is  possible  for  (HI)  to  give  a schedule  which  not  only  Is  not 
optimal,  but  which  can  be  easily  Improved.  This  section  characterizes  such 
situations  and  proposes  a collapsing  routine  which  Is  to  be  applied  to  the 
schedule  generated  by  (HI) . 

Consider  a two-stage  system  (l.e.  Figure  2),  and  assume  there  is  no 
external  demand  for  item  2 (d 2t  “ 0) . For  this  system  it  is  quite  possible 
for  (HI)  to  give  a schedule  such  that 

(5)  x"t  - ex"t  for  t-1,2 T. 

That  is.  Item  2 Is  simultaneously  produced  whenever  Item  1 is  produced,  and 
no  Inventory  Is  ever  kept  for  item  2.  Note  that  If  (x”t)  is  the  schedule 
for  item  1 after  the  nth  Iteration  of  (HI),  and  If  x^  ” 6x"t>  then  we 
have  x"*1  “ x”t  which  implies  convergence  for  (HI) . This  will  occur 
because  at  iteration  n+1  the  "revised"  variable  production  cost  of  Item  1 
In  product  t Is  unchanged  If  x”fc  > 0,  but  is  higher  If  x”£  ” 0.  Consequently, 
the  new  schedule  for  Item  2 does  not  alter  the  scheduling  of  item  1. 

An  implication  from  (5)  is  that  the  two-stage  problem  may  be  restated 
as  a single-stage  WW  problem  with  setup  cost  S - + S^,  holding  cost 

h ■ h^,  variable  production  cost  v£  ■ v^fc  + 6v2t»  and  demand  dfc  ■ dlfc. 

The  otpimal  solution  of  this  single-stage  problem  must,  by  definition,  be 
no  worse  than  the  schedule  in  (5).  Thus,  if  (HI)  generates  a schedule 
satisfying  (5),  by  collapsing  stage  2 into  stage  1 we  can  define  a single- 
stage  problem,  from  which  an  improved  schedule  may  be  found. 

This  improvement  can  be  generalized  for  more  complex  N-stage  problems. 
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We  need  restate  condition  (5)  for  stage  1 as 

(6)  m t“l,2, . . . ,T, 

where  stage  J Is  the  unique  successor  to  stage  1.  Now  stage  1 may  be 
collapsed  Into  stage  J by  redefining  stage  j to  have  setup  cost  s “ Sj  + s^> 

holding  cost  h - hj , and  variable  production  cost  v£  - Vj£  + BjjV^.  Note 
that  the  immediate  predecessors  or  Inputs  to  stage  1 are  now  direct  Inputs 
to  stage  j.  After  collapsing  as  many  stages  as  possible,  we  may  now  reapply 
heuristic  (HI)  to  find  an  Improved  schedule;  the  procedure  stops  when  a 
schedule  is  found  which  cannot  be  Improved  by  (HI) , and  from  which  no  more 
stages  can  be  collapsed. 


3.  Computation  Tests 


In  order  to  test  the  performance  of  the  proposed  multipass  heuristic, 

we  tested  the  heuristic  procedure  on  50  test  problems  generated  for  each  of 

five  multistage  assembly  systems.  Figure  5 gives  the  five  assembly  systems 

each  with  five  stages;  these  systems  are  the  same  as  used  in  [11],  [12].* 

In  generating  the  test  problems,  we  normalized  h^  ■ 1;  for  j-l,...,4  we  set 

h - e + l h where  B(j)  is  the  set  of  immediate  successors  to  j and 
J i£B(j) 

e is  selected  from  a uniform  distribution  with  values  e » 0.1,  0.5,  1.0, 
and  2.0.  We  selected  from  a uniform  distribution  with  values  « 150, 
300,  600,  1500;  we  selected  dlt  from  a uniform  distribution  with  values 
dlt  - 0,  10,  20,  30,  40,  100,  200,  400.  There  was  no  other  external  demand 
(d^t  “ 0 for  all  t,  and  i-2,...,5),  all  variable  production  costs  v^t  were 
zero,  and  the  length  of  the  horizon  T was  12  periods. 

For  each  test  problem,  we  compared  the  multipass  heuristic  with  two 
single-pass  heuristics  and  with  an  optimal  algorithm.  The  mutlipass 
heuristic  (MP)  that  was  implemented  for  these  test  problems  iterated  in  a 
full  cyclic  fashion  over  the  stages;  that  is,  the  procedure  scheduled 
stages  1,2,..., 5 using  a Wagner-Whltln  algorithm,  and  then  returned  to 
stage  1 to  begin  the  revisions.  Alternative  schemes  were  examined,  but  did 
not  seem  to  give  significant  improvements.  The  optimal  schedule  was  found 
using  the  dynamic  programing  procedure  of  Crowston  and  Wagner  [4].  One 
single-pass  heuristic  (SP-WW)  was  to  schedule  each  stage  using  the  Wagner- 
Whltln  algorithm.  The  second  single-pass  heuristic  (SP-MW)  was  that  proposed 

* Although  all  the  test  problems  are  for  assembly  systems,  the  multipass 
heuristic  can  be  used  for  any  multistage  system.  Assembly  systems  were 
chosen  for  the  computational  tests,  because  the  multipass  heuristic  is  to 
be  compared  with  a single-pass  heuristic  proposed  in  [12]  and  with  an  opti- 
mal algorithm  from  [4],  both  of  which  are  restricted  to  assembly  systems. 
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Figure  5:  Multistage  Assembly  Systems 
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by  McLaren  and  Whybark  (121;  here  each  stage  Is  scheduled  using  the  Wagner- 
Whitin  algorithm,  but  with  the  stage's  setup  cost  inflated  to  reflect  the 
possible  setups  for  preceding  stages. 

The  primary  results  are  given  in  Tables  1 and  2;  Table  1 reports  the 
number  of  times  each  heuristic  obtained  the  optimal  solution,  while  Table  2 
gives  the  average  percentage  cost  penalty  for  the  heuristics  over  the 
optimal  procedure.  From  these  tables,  in  the  aggregate  the  SP-WW  is  domin- 
ated by  the  SP-MW  which  is  dominated  by  the  MP.  However,  for  each  problem 
structure  there  were  exceptions  to  this  domination  such  that  on  a few  test 
problems  SP-WW  outperformed  SP-MW,  or  SP-MW  outperformed  MP;  SP-WW  never 
did  better  than  MP,  since  the  multipass  heuristic  is  an  improvement 
routine  which  is  initiated  with  the  SP-WW  schedule.  We  should  note  that 
even  though  the  MP  is  optimal  in  about  90Z  of  the  test  problems,  compared 
with  a rate  of  about  60Z  for  SP-MW,  both  heuristics  have  minuscule  average 
percentage  cost  errors.  Also,  the  performance  of  both  heuristics  seemed 
to  vary  with  the  number  of  echelons  in  the  system;  the  performance  relative 
to  the  optimal  procedure  becomes  poorer  as  the  number  of  echelons  grows. 

The  cooputational  effort  required  for  the  heuristics  and  for  the 
optimal  procedure  was  consistent  with  their  respective  performances.  Both 
the  SP-WW  and  the  SP-MW  heuristics  solve  one  WW  problem  per  stage,  or  five 
WW  problems  for  each  test  problem;  on  a PRIME  400  minicomputer,  the  250 
test  problems  using  the  SP-MW  heuristic  were  solved  in  53  seconds  of  cpu 
time  (0.21  seconds/problem).  The  multipass  heuristic  due  to  its  iterative 
nature,  must  solve  a variable  number  of  WW  problems.  It  can  be  shown  that 
a gross  upper  bound  on  the  number  of  WW  problems  considered  by  the  multipass 
heuristic,  as  Implemented  in  this  work,  is  T*N  where  T is  the  length  of  the 
horizon  and  N is  the  number  of  stages;  since  the  WW  problem  can  be  solved 
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Table  1:  Frequency  of  Optimal  Schedules  (out  of  50  problems) 

HEURISTIC 


SP-WW 

SP-MW 

MP 

A 

2 

19 

41 

B 

4 

42 

48 

C 

2 

27 

44 

D 

2 

20 

46 

E 

4 

31 

47 

2:  Average  Percentage  Cost  Errors 

HEURISTIC 


SP-WW 

SP-MW 

MP 

A 

8.64Z 

2.10Z 

0.58Z 

B 

3.27Z 

0.06Z 

0.05Z 

C 

5.12Z 

0.64Z 

0.26Z 

D 

7.45Z 

1.70Z 

0.31Z 

E 

4.76Z 

0.61Z 

0.26Z 

efficiently,  the  multipass  heuristic  is  an  efficient  procedure.  For  the 
250  test  problems,  the  average  number  of  WW  problems  considered  was 
slightly  over  18  per  test  problem;  the  total  cpu  time  for  the  250  test 
problems  was  115  seconds  (0.46  seconds/problem).  The  optimal  algorithm 
as  reported  in  [4  ],  is  an  inefficient  procedure,  in  that  the  computational 
work  grows  exponentially  with  the  length  of  the  horizon;  for  the  250  test 
problems  the  optimal  algorithm  required  2077  seconds  of  cpu  time  (8.31 
seconds/problem) . 

In  summary,  the  multipass  heuristic  performs  better  than  the  leading 
single-pass  heuristic,  SP-MW,  but  requires  more  work.  In  light  of 
the  computational  experience,  we  conjecture  that  the  effort 
needed  for  the  multipass  heuristic  will  be  proportional  to  that  needed 
for  a single-pass  heuristic.  Finally,  we  note  chat  whereas  the 
multipass  heuristic  is  applicable  to  general  multistage  systems,  both  the 
SP-MW  heuristic  and  the  optimal  algorithm  are  designed  for  assembly  systems. 
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APPENDIX:  Proof  of  Convergence  for  Three-Stage  Serial  System 

Consider  the  three-stage  system  depicted  In  Figure  4.  Assume  that  the 
revision  procedure  as  outlined  in  Section  2.2  has  been  applied  where  the 
iteration  over  the  echelons  is  in  a pure  cyclic  fashion  (e.g.  echelon 
1 ■*  II  ■*  III  -*•  I -*•  II  . , . ) . To  show  that  the  procedure  converges,  we  will 
follow  the  arguments  given  in  Section  2.1  for  the  two-stage  system.  For 


notational  ease  we  assume  that  0 
n+1  -i 


32 


$22  * 1.  After  the  n+lSt  iteration, 


for  (x,.  } being  the  schedule  for  Item  1 with  actual  cost  c,  the 

“ T l,n+l 

"revised"  cost  using  v is  given  by  c.  . + £ Y^x!**1.  Since  {x"*1} 

it  1 iIItI  . it  It  It 

t-i 

is  optimal  for  the  "revised"  problem,  we  have 

T T 

/ , , \ I ^ n n+ 1 . . « n n 

(A1)  cl,n+l  + Z Y2txlt  i cln  + Z.  Y2tXlt 

t“l  t“l 

or 


(A2) 


where  Ax,_ 

it 


Cl,n+1  1 Cln  ' tZml  Y2t(Axlt) 


n+1  n 
Xlt  " xlt' 


For  item  2 at  the  n+lSt  iteration,  the  demand  is  given  by  “d2t+  Xl’t^‘ 

A feasible  schedule  { y t } may  be  derived  from  the  previous  schedule  {x^}  as 
follows: 


and 


2t 


if 

n 

2t 

T— 1 

l 

Ax, , 

1J 

0, 


if  x“t  > 0 


where  T is  defined  as  the  earliest  period  after  t such  that  x^.  > 0.  The 
cost  associated  with  {y£}  for  the  "revised"  problem  ^l^t '^2 ’^2 ,V2t ^ *8 


(A3) 


c2o + ^ 


\ 


- . -■sVV,*.. 


- i irtfci  dw 
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where  c„  is  the  actual  cost  associated  with  schedule  (x!)  }.  To  verify 
2n  T 2t 

(A3),  note  that  c_  + E y”  x"  is  the  "revised"  cost  for  schedule  (x"  } 

2n  3t  2t  ^t 

assuming  demand  d"t  “ d£^  + xlt‘  ^en  the  item's  demand  is  restated  as 
d"+*  “ d!J  + Ax.  , schedule  {x"  } need  not  be  feasible;  the  feasible  schedule 

it  It  Zt  j 

{y  } dervied  from  {x”  },  has  additional  costs  which  are  given  by  E y_  (Ax^  ) , 
t -t  4t 

by  definition  of  the  marginal  costs. 

The  cost  for  the  optimal  schedule  for  the  "revised"  problem 

WW[d2t.S2*h2»''2t^  ls 


(A4) 


v n n+1 

C2,n+1  Z 73tX2t 
t*l 


where  c.  ,,  is  the  actual  cost  of  the  schedule.  Due  to  the  optimality 
L , n+1 

of  (x^>  we.  have  that  (A4)  is  less  than  or  equal  to  (A3),  which  may  be 
written  as 


(A3) 


C2,n+1  i C2n  + ^ 4(Axlt)  " ^ '"t'^t^ 
s t 

For  item  3 at  the  n+1  iteration,  the  demand  is  given  by 


d”*1  - djt  + x2^‘  T^e  Previous  schedule  {x"t}  may  be  modified  to  satisfy 
this  demand  in  a similar  fashion  to  the  derivation  of  {ytJ  for  item  2.  The 


actual  cost  for  this  schedule  can  be  shown  to  be  equal  to  c^n  + I , >3t(Ax2t> 

n+li  C 1 

C3t 


If  schedule  {x”  } with  cost  c,  . is  the  optimal  schedule  to 


‘3, n+1 

WWld^.Sj.hj.v^] , then  we  must  have 

(A6)  c3,n+l  ± C3n  + ^ Y5t(Ax2t) 

due  to  the  optimality  of  {x"^}« 

By  combining  the  inequalities  (A2) , (A5) , and  (A6) , we  have 


(*  + p + r < c 

l,n+l  2, n+1  3,n+l  - In 


"2n 


'3n* 
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which  states  that  at  each  iteration  total  cost  is  nonincreasing.  By 
similar  reasoning  to  that  given  for  the  two-stage  example,  this  result  is 
sufficient  to  guarantee  convergence. 

To  extend  this  result  to  an  N-stage  serial  system  is  straightforward. 
Inequalities  (A2)  and  (A6)  can  be  shown  true  for  stage  (item)  1 and  for 

stage  (item)  N,  respectively.  For  stages  2,3 N-l,  an  inequality 

identical  to  (A5)  can  be  established.  By  combining  these  inequalities  the 
desired  result  is  obtained. 

For  more  general  multistage  systems,  the  convergence  proof  is  similar 
in  structure  but  is  complicated  by  the  notational  needs.  Essentially,  for 
each  echelon  in  the  system  an  inequality  similar  to  (A2),  (A5),  or  (A6) 
is  derived  which  relates  the  total  actual  costs  of  the  echelon  at  the 

S t t" 

n+1  iteration  to  those  costs  at  the  n iteration.  The  details  of  this 
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